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Identifying the sites of r-process nucleosynthesis, a primary mechanism 

of heavy element production, is a key goal of astrophysics. The discovery 
of the brightest gamma-ray burst (GRB) to date, GRB 2210094, presented 
anopportunity to spectroscopically test the idea that r-process elements 
are produced following the collapse of rapidly rotating massive stars. Here 
we present James Webb Space Telescope observations of GRB 221009A 
obtained +168 and +170 rest-frame days after the gamma-ray trigger, and 
demonstrate that they are well described by a SN 1998bw-like supernova 
(SN) and power-law afterglow, with no evidence for a component from 
r-process emission. The SN, with a nickel mass of approximately 0.09 M,, is 
only slightly fainter than the brightness of SN 1998bw at this phase, which 
indicates that the SN is not an unusual GRB-SN. This demonstrates that the 
GRB and SN mechanisms are decoupled and that highly energetic GRBs 

are not likely to produce significant quantities of r-process material, which 
leaves open the question of whether explosions of massive stars are key 
sources of r-process elements. Moreover, the host galaxy of GRB 221009A 
has a very low metallicity of approximately 0.12 Z, and strong H, emission 
at the explosion site, which is consistent with recent star formation, hinting 
that environmental factors are responsible for its extreme energetics. 


The origin of the heaviest elements in the Universe, specifically those 
formed by means of rapid neutron capture (r-process) nucleosynthe- 
sis, remains a major open question in astrophysics!”. Given the high 
density of neutron-rich material needed for the r-process to occur, the 
collisions of neutron stars have long been a suspected source™* and, 
indeed, the observations of the kilonova associated with GW 170817 
confirmed that binary neutron star (BNS) mergers are the source of 
atleast some of the r-process material in the Universe^?. However, 


thereisgrowing evidence that there may be multiple sites of r-process 
nucleosynthesis from studies of low-metallicity galactic halo stars, 
dwarf galaxy and globular cluster enrichment". 

Asecond proposed site of the r-process is in rapidly rotating cores 
of massive stars that collapse into an accreting black hole, producing 
similar conditions as the aftermath of a BNS merger”. Theoretical 
simulations suggest that accretion disk outflows in these so-called 
‘collapsars’ may reach the neutron-rich state required for the r-process 
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Fig. 1]JWST/NIRCam imaging. Our images of GRB 221009A (top row), best-fit 
GALFIT galaxy models (middle row) and GALFIT model subtracted images 

(bottom row). Images are shown with north up and east to the left. A clear point 
source is detected at the location of GRB 2210094. The red rectangle shows the 


NIRSpec slit orientation. PSF photometry of GRB 221009A was performed on 
the galaxy-subtracted images. The host galaxy is well described by a single Sérsic 
component, although some residual galaxy structure remains in the F200W, 
F277W and F444W filters. 


to occur ^*^, The larger mass of r-process material synthesized per 
event compared with BNS mergers suggests that collapsars could 
be a dominant source, making them a possible missing piece in our 
understanding of r-process enrichment in the Universe. 

The discovery ofthe long-duration gamma-ray burst (GRB) GRB 
221009A, the brightest GRB ever observed" ?, on 9 October 2022 at a 
relatively nearby redshift of z= 0.151 (ref. 20) presents a unique oppor- 
tunity to search for r-process signatures in a collapsar. Collapsars are 
the favoured explanation for long GRBs (LGRBs), which result from 
the launch of a relativistic jet and its subsequent interaction with the 
surrounding medium"? r-Process nucleosynthesis is more likely to 
occur in collapsars with large accretion disk masses, which are also 
thought to be linked with brighter GRBs”, making GRB 221009A a 
particularly strong candidate to search for r-process signatures. These 
events are known to be accompanied by broad-lined type Ic supernovae 
(SNe Ic-BL) characterized by higher velocities than normal type Ic 
supernovae, suggesting that the energy powering LGRBs also affects 
the associated supernovae (see ref. 24 for a review). 

It is the supernova (SN) following aLGRB that would be responsible 
for carrying r-process material from the explosion site into the inter- 
stellar medium. Although early-time observations of GRB 221009A 
provided an exquisite view of the afterglow” 7, to date, there are con- 
flicting claims in the literature regarding the presence of an associ- 
ated SN, which are due, in part, to the bright afterglow and high Milky 
Way extinction? ?', Moreover, there have been claims that two recent 
LGRBs are associated with BNS mergers? ?, making the search for 
an SN associated with GRB 221009A crucial not only for an r-process 
search, but also for understanding the origin of its extreme luminosity. 

Here, we present late-time James Webb Space Telescope (JWST) 
observations of GRB 2210094 consisting of a near-infrared (NIR) spec- 
trum and imaging in four NIR bands. These observations provide clear 
detection of an SN associated with this extreme event and enable the 
search for r-process emission in a nebular-phase spectrum of a GRB-SN. 
Moreover, these data provide a detailed NIR view of the host galaxy, 
enabling an assessment of environmental factors that may be respon- 
sible for this extraordinary GRB. 


Identification of SN emission 
We obtained spectroscopy with the Near Infrared Spectrograph 
(NIRSpec) using the medium-resolution gratings covering 1-3 um on 


20 April 2023 and imaging with the Near Infrared Camera (NIRCam) 
using the F115W, F200W, F277W and F444W filters on 22 April 2023. 
These observations occurred +194 and +196 observer-frame days 
after the burst (rest-frame phases of +168 and +170 days, respectively). 
The afterglow of GRB 2210094 is clearly detected in our images, from 
which we measured photometry (Fig. 1; see Methods for details). In our 
NIRSpec observations, we detecta clear spectral trace containing flux 
from GRB 221009A and its host galaxy (Extended Data Figs. 1and 2; see 
Methods for details ofthe spectral extraction). Owing to the high Milky 
Way extinction and possible non-negligible extinction intrinsic to the 
host galaxy??? we analysed archival early-phase NIRSpec/PRISM and 
Mid-Infrared Instrument (MIRI) spectra? of GRB 2210094 using multi- 
ple dust laws to constrain the extinction (Methods and Extended Data 
Fig. 3). We used the resulting extinction parameters (Extended Data 
Table 1) to correct our rest frame +168 day NIRSpec grating spectrum. 

In Fig. 2, we show two versions of the spectrum, one corrected 
using an extinction curve from ref. 37 and another using one from ref. 
38, transformed to the rest frame of GRB 221009A. In both cases, the 
spectrum exhibits an overall flat shape in the range -1-1.5 um, witha 
smooth, gradual upturn at redder wavelengths extending to the edge 
of our coverage at -2.7 jum and a sharp upturn at bluer wavelengths 
due, in part, to apparent broad emission features. The use of different 
extinction laws and parameters, within the range of uncertainties from 
our fitting, does not change these fundamental characteristics. 

The gradual rise in the spectrum at wavelengths 2 1.5 um strongly 
resembles a power-law shape and therefore this region is likely to 
contain a significant contribution from the afterglow of GRB 221009A. 
In addition, our photometric observation in the F444W filter (which 
lies redward of our spectroscopic coverage) indicates that the flux 
continues to rise at longer wavelengths (23.8 um, rest frame). The fluxes 
measured in the F200W, F277W and F444W filters are consistent with 
a single power law with an index of f = -0.64 + 0.10. The shape of the 
NIRSpec spectrum at wavelengths A 21.5 um is slightly steeper than 
this slope, with a power-law index of B = -0.76 + 0.07, although it is 
consistent within the uncertainties. 

AtA € 1.5 um, the spectrum clearly deviates from an extrapola- 
tion of the power law at A z 1.5 um, exhibiting an overall flat shape and 
several broad SN-like emission features. Indeed, we identify two broad 
emission features located at wavelengths of -0.86 jum and -0.92 um, 
which are consistent with the Ca II NIR triplet and O I, respectively. 


Nature Astronomy 


Article 


https://doi.org/10.1038/s41550-024-02237-4 


4.0 
—— NIRSpec corrected with G23 law (A, = 4.37, R, = 3.07) 
35 4 —— NIRSpec corrected with F99 law (A, = 4.64, R, = 4.24) 
€ NIRCam corrected with F99 law 
304 —— ALMA-XRT power law at 168 rest-frame days 


Flux density (uJy) 


Rest wavelength (um) 


Fig. 2 |JWST/NIRSpec spectrum. Our +168 rest-frame phaseJWST/NIRSpec 
G140M+G235M spectrum of GRB 221009A corrected for extinction (see Methods 
for spectral extraction details; error spectra are shown in Extended Data Figs. 1 
and 2). We show two versions corrected using the G23 (orange)? and F99 (blue) 
extinction laws and corresponding best-fit extinction parameters from fitting 
the early-time NIRSpec/PRISM and MIRI data from ref. 31 as described in Methods 
"Constraints on foreground dust from early-time spectroscopy’. In both cases, 
thespectrum appearsto exhibit multiple components, with SN-like emission 
atA<1.5 pmand rising flux at A 2 1.5 jum, which is likely to be due to the GRB 
afterglow power law. We also show ourJWST/NIRCam photometry corrected 
using the F99 extinction law” (points with lo error bars), as well as the ALMA-XRT 
power law (black line). 


These are typical nebular-phase emission lines observed in core- 
collapse supernovae (for example, ref. 39). We show a zoomed-in 
comparison of these features with those seen in SN 1998bw, SN 
2013ge and SN 2014ad in Extended Data Fig. 4. In addition to the flat 
spectral shapeat -1-1.5 um, these emission features strongly support 
the identification of SN emission in our spectrum of GRB 2210094. Our 
observation, therefore, represents the latest phase NIR spectrum of an 
SNassociated with a GRB to date. 


Isolating the SN signal 

Although disentangling the SN and afterglow components is 
not straightforward, the relative featureless nature of the red end of 
the spectrum indicates that the afterglow component is sufficiently 
bright to not only affect the overall shape but also to dilute SN 
features with respect to the continuum in that region (see Methods 
and Extended Data Fig. 5 for comparisons with previous supernovae). 

To separate the afterglow and SN components, we considered 
several afterglow models. First, we used Atacama Large Millimeter/ 
submillimeter Array (ALMA) and Swift X-ray Telescope (Swift/XRT) 
observations obtained at roughly the same phase as our NIRSpec spec- 
trum and modelled the afterglow at NIR wavelengths as a power law 
connecting the radio and X-ray data. We find flux density F, « y °° *°% 
for frequency v (see Methods for details). We show this power law, 
normalized to the measured radio and X-ray flux, compared with our 
spectrum in Fig. 2. While the ALMA-XRT power-law slope is similar 
to the shape of our spectrum at A 21.5 uum, our data is systematically 
offset to higher flux, which indicates that the ALMA-XRT power law 
does not fully capture the afterglow contribution at NIR wavelengths. 
Moreover, the implied SN component deviates from the expected 
spectral shape of an SN (Methods and Extended Data Fig. 5). 

Next, we modelled the afterglow from our spectrum itself, namely, 
as a power law with a slope determined by fitting our spectrum at 
wavelengths z 1.5 um, where the afterglow is likely to be dominating. 
We find a best-fit power law of F, « y ?76*09?7 This is steeper than the 
ALMA-XRT power law, which further confirms that interpolating the 
millimetre and X-ray bands is not likely to provide the best representa- 
tion of the afterglow at these wavelengths. We then performed a joint 
fit ofan SN template and the fitted power law, with the power-law slope 


d € 401 PES — SN 1998bw «51/473 scaled to +168 d 
3 304 pared ok 5s —— SN 2013ge +118 d scaled 
i 204 — Best-fit SN + afterglow | — sn 2014ad +79 d scaled 
10 d a 
E 2 
j © 104 
a ke] d 
H x i 
=i 
T [ras 


1.0 1.5 2.0 25 


10° 
Rest wavelength (um) 


Flux density (uJy) 


i i it ' i | 


| m 


2.00 25 


e TW nm AMA TN 1 
IP l Midh 
Afterglow model | Ti ll 


subtracted 
T 


T T 
0.75 1.00 1.25 1.50 1.75 2.50 2.75 


Rest wavelength (um) 


Fig. 3 | Spectral fit and comparisons. Our NIRSpec spectrum of GRB 

221009A (blue, smoothed) after subtracting our best-fit afterglow model. The 
unsubtracted spectrum, best-fit afterglow model (dashed black line) and best-fit 
SN + afterglow (black line) are shown in the inset. We show late-time spectra of SN 
2013ge (gold?) and SN 2014ad (magenta?) scaled to match the shape and features 
oftheafterglow-subtracted spectrum atA € 1.5 im where the SN dominates, 
demonstrating the overall resemblance with these comparison supernovae Ic/ 
Ic-BL. We also show SN 1998bw (orange*?) scaled to the distance of GRB 221009A 
andthe phase of our spectrum, showing that it matches not only the shape but 
the overall flux level of our spectrum. The close match with supernovae Ic-BL, in 
particular, demonstrates the presence of a typical GRB-SN in our spectrum. 


fixed, to determine the best-fit combination of SN and afterglow. For 
the SN template, we used the +51 day spectrum of SN1998bwas this is 
the latest available NIR spectrum of another GRB-SN^?, allowing the 
overall flux normalization to vary. 

In Fig. 3, we present the best-fit SN 1998bw + afterglow spectrum 
and our spectrum of GRB 221009A after subtracting the best-fit after- 
glow component. We compare our afterglow-subtracted spectrum 
with the SN 1998bw spectrum scaled to the distance of GRB 221009A 
and the brightness of SN 1998bw at the phase of our JWST spectrum 
using the light curve of SN 1998bw from ref. 41. The best-fit SN com- 
ponent is ~30% fainter than the expected brightness SN 1998bw would 
have at this distance and phase. We also compare with late-time spec- 
tra of the SN Ic SN 2013ge? and the SN Ic-BL SN 2014ad”. To directly 
compare the shapes andfeatures we scaled SN 2013ge and SN 2014ad 
to best match the spectral shape and features at the blue end of the 
afterglow-subtracted spectrum where the SN component dominates. 

These events provide an excellent visual match to the afterglow- 
subtractedspectrum, which confirms that our estimate ofthe afterglow 
contributionis reasonable. In addition, the inferred ratio of Ca 11/0 ris 
amuchbetter matchto theratios seeninthethree comparison objects 
compared with the case of no afterglow subtraction (Extended Data 
Fig. 4). Although the width of the Ca 11 emission complex exhibits a 
better match with SN 2013ge, the afterglow-subtracted spectrum does 
not show the same strong absorption seen at -1.1 uum in SN 2013ge, 
possibly due to the SN associated with GRB 2210094 having a higher 
ejecta velocity. SN Ic-BL-like velocities are further supported by the 
better overall match to SN 2014ad and SN 1998bw. The narrower width 
of Ca ri compared with SN 2014ad and SN 1998bw may be an artefact 
of the instrumental response affecting the shape at the blue end of 
the line. We also identify evidence for a broad emission feature near 
A=1.5 um, which is consistent with the location of the 1.503 um line 
of MgIseen in the comparison objects and large samples of other 
supernovae Ic/Ic-BL?. 

In summary, our spectrum is well fit by an SN and power-law 
model; we do not require another component to explain the spectrum, 
although we explore the possibility that the afterglow contribution is 
lower and whether some of the resulting red excess in such a model 
(Methods section ‘Constraints on the afterglow contribution’) could 
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be explained by r-process emission in the section ‘No sign of r-process 
enrichment’ and Methods section ‘Comparison with r-process light 
curve models’. Importantly, our afterglow-subtracted spectrum is 
similar to, although slightly fainter than, the expected flux of SN1998bw 
at the distance of GRB 221009A and the phase of our observations, 
which suggests that the SN associated with GRB 221009A produced a 
similar quantity of *°Ni. 


A modest nickel mass indicates a typical GRB-SN 

Estimates for the mass of *°Ni produced in SN 1998bw range from 
-0.3 M, to -0.7 M, depending on the models and assumptions 
used to fit the light curve“. Reference 47 considered a two-zone 
model where «0.44 M, of?9Ni is contained in an outer high-velocity 
componentthat rapidly expands and becomes optically thin, explaining 
the bright peak luminosity. An additional «0.12 M, exists in an inner 
dense low-velocity component that explains the linear nature of the 
light curve at intermediate phases of -100-200 days. 

We directly estimated the mass of *°Ni produced by the SN asso- 
ciated with GRB 2210094 by integrating the afterglow-subtracted 
spectrum. We estimated the unobserved flux using SN 2007gr as a 
spectral template owing to its simultaneous optical and NIR coverage 
outtothe same phase of our observations. We find thatthe wavelength 
coverage of our NIRSpec spectrum accounts for about 50% of the 
total emitted flux. At the phase of our observations, the luminosity of 
a nickel-powered SN is dominated by the decay of °Co, the daughter 
isotope of Ni. Assuming a single component of the ejecta and full 
gamma-ray trapping, we find My; ~= 0.03 M,. Under a more realistic 
assumption of gamma-ray leakage, with a timescale of -100 days for 
the ejecta to become optically thin to gamma-rays (as inferred for 
SN1998bw), we find My; = 0.09 M. 

The *°Ni mass we infer assuming gamma-ray leakage is therefore 
slightly lower than the mass inferred by ref. 47 for the inner dense 
component of SN 1998bw, which is consistent with our inference 
that the SN associated with GRB 2210094 is slightly fainter than SN 
1998bw at late time. Of course, assuming a different afterglow contri- 
bution in our spectrum will affect the estimated mass. Our inferred 
mass is consistent with the results of ref. 30, who foundbest-fit values 
from modelling the early light curve of GRB 2210094 in the range 
My; = 0.05-0.25 M,, depending on assumptions about the host extinc- 
tion, with a 99% upper limit of My; < 0.36 M,. These values are lower than 
most early light curve estimates for SN1998bw***>”, This may indicate 
alower ratio of the outer-to-inner ejecta components compared with 
SN 1998bw, or that a two-component model is not needed to explain 
the SN associated with GRB 2210094. Our results, combined with the 
early light curve estimates, conclusively rule out the possibility that 
theSN was unusually bright compared with previous GRB-SNe. This is 
consistent with previous sample studies that do not show a correlation 
between the luminosities of LGRBs and their associated supernovae^*^, 
Crucially, our spectroscopic detection of the SN confirms that the 
marginal deviation from a typical afterglow in the early light curve 
claimed by ref. 29 and ref. 30 was indeed due to the SN. 


Nosigns of r-process enrichment 

The identification of the SN associated with GRB 2210094 allows us 
to constrain the presence of r-process material. One possibility is that 
the red excess in our spectrum consists of acombination of afterglow 
and emission from r-process elements. Reference 15 outlines how 
acollapsar with a massive transient disk can lead to r-process produc- 
tion. However, the observational impact of r-process material, if it 
is produced, is highly dependent on the degree of outward mixing. 
In particular, ref. 15 presents two possible scenarios: one (the 'magneto- 
hydrodynamic (MHD)' case) in which 0.025 M, of r-process elements are 
mixed uniformly throughout the SN ejecta with v « 0.15c and one (the 
'collapsar' case) in which 0.25 M, of r-process elements are confined 
tov « 0.015c. In both cases, 0.25 M, of?9Ni is mixed in the ejecta. In the 
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Fig. 4 | Comparison with r-process models. Comparison of our NIRSpec 
spectrum of GRB 221009A with r-process enriched SN models from ref. 15 
corresponding toa phase of 95 days after explosion, which is the latest phase 
available. We show our original spectrum without afterglow subtraction (top 
blue), as well as the resulting spectra after subtracting the ALMA-XRT power law 
(middle black) and our best-fit afterglow model (shown in Fig. 3, bottom, light 
blue). Our spectrum, even after accounting for the afterglow, is clearly distinct 
from the predictions of an MHD SN. We also do not see evidence for spectral 
features in our spectrum that can be linked to the collapsar wind model and not 
attributed to the SN. 


MHD case, the r-process material tracks 5Ni, whereas in the collapsar 
case, the r-process elements are embedded behind the **Ni. In truth, 
the degree of mixing in the collapsar wind scenario is unknown, andis 
likely to be variable with progenitor properties and may be sufficient 
to mix r-process elements with the outer layers. 

Although the MHD scenario has largely been ruled out by 
early-time observations of previous events’, few constraints exist 
onthe collapsar wind scenario dueto the lack of late-time NIR spectra 
of GRB supernovae. Before our NIRSpec spectrum of GRB 221009A, 
the latest NIR spectrum of a GRB-SN was that of SN 1998bw taken at 
+51 days, which we have shown is an excellent match to our spectrum 
(Fig. 3) after subtracting our best-fit afterglow power law. Here we 
consider the possibility that our best-fit power law overestimates the 
afterglow contribution and that our much later spectrum of the SN 
associated with GRB 221009A differs from the +51 day NIR spectrum 
of SN 1998bw owing to the presence of r-process signatures. 

In Fig. 4, we compare our NIRSpec spectrum, with various assump- 
tions about the afterglow contribution, to r-process enriched SN 
models from ref. 15 (with r-process masses up to 0.25 M,). We compare 
with models corresponding to a phase of 95 days after explosion, the 
latest phase available, and shift them to the distance of GRB 221009A. 
At this phase, the MHD SN differs considerably from an SN without 
r-process enrichment, producing strong emission at «1.8-2.4 um that 
is clearly not present in our spectrum whatever the assumption on 
afterglow contribution. The collapsar wind model, onthe other hand, 
largely shows SN features from non-r-process elements, although 
with enhanced flux near «2 um compared with what is seen in normal 
supernovae. 

Owing to the noise in our spectrum, we are unable to identify 
individual lines in this region of the spectrum, beyond the likely Mg I 
atA«1.5 um. However, we compared the overall flux level and find that 
assuming no afterglow contribution (thatis, the original unsubtracted 
spectrum) leads to much higher continuum flux than the collapsar wind 
model forA > 1 um. Furthermore, the expected strong nebular SN lines 
are diluted with respect to the continuum (Methods), indicating an 
extra continuum source is present (the afterglow of the GRB). Assuming 
the afterglow shape and normalization given by interpolating the 
contemporaneous ALMA and XRT observations, we find that, overall, 
the spectrum is inconsistent with the collapsar wind model, which 
indicates that the spectrum cannot be explained by a combination 
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Fig. 5 | Host galaxy spectra. Spectrum of the total host galaxy including the 
site of the GRB (grey), the "host-only' spectrum excluding the GRB site (red; 
see Methods section ‘Late-time NIRSpec observations’) anda narrow aperture 
centred on the location of GRB 2210094 (blue). We detect narrow H, H;, He, 
FeandS emission lines from the galaxy. Importantly, we see that some narrow 


emission lines change in strength over this galaxy, having notably strong 
molecular H, emission in the region of GRB 221009A. The continuum ofthe 
afterglow and SN component canbe seen clearly in the blue spectrum asa 
deviation from the host-only spectrum. 
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Fig. 6 | Best-fit galaxy model. Left, best-fit Prospector model photometry (green circles) and spectrum (green) compared with the observed photometry shown with 
lcerror bars (orange circles, our NIRCam photometry and HST photometry measured by ref. 31). Right, best-fit Prospector spectrum compared with our NIRSpec 


spectrum (orange). The bottom panels present the residuals. 


of the ALMA-XRT power law and an r-process enriched SN. Given the 
strong resemblance to previous supernovae across the full wavelength 
range when assuming our best-fit afterglow power-law shape and 
contribution (Fig. 3), itis unlikely that flux from r-process elements are 
contributing significantly to our spectrum. Our observation highlights 
the need for a systematic survey of nebular-phase LGRB spectra across 
a broad range of GRB properties, in particular, in light of the recent 
theoretical work that correlates these properties to the degree of 
r-process production. We additionally compare our observations 
to the broadband colour evolution models due to r-process enrich- 
mentfrom ref. 52, which further highlights the need for spectroscopy 
(Methods section 'Comparison with r-process light curve models"). 


Host galaxy properties 

A very low-metallicity, star-forming galaxy 

The host galaxy of GRB 221009A is readily apparent in our JWST/ 
NIRCam imaging shown in Fig. 1. Consistent with analysis of the opti- 
cal Hubble Space Telescope (HST) images”, we find that GRB 221009A 
is located 0.24 + 0.01 arcsec (0.66 + 0.02 kpc assuming the cosmo- 
logical parameters presented in ref. 53) from the centre of its host 
galaxy, which appears to be a near edge-on system. From our GALFIT 


modelling (Methods) we find that this galaxy is well described by a 
single Sérsic component with index n = 1.2 + 0.1 and effective radius 
r, = 2.15 + 0.07 kpc. These values represent the mean and standard 
deviation across the four filters. The AB magnitudes in each filter cor- 
responding to the best-fit GALFIT models are mg, = 21.58 + 0.20 mag, 
Ilig;9ow = 20.62 + 0.10 mag, Mp277w = 20.88 + 0.10 mag and Mrp444w = 
21.38 + 0.05 mag (not corrected for Milky Way extinction). 

In Fig. 5, we show the global host spectrum (that is, including 
flux from the entire resolved spectral trace) and the spectrum at the 
position of GRB 221009A (see Methods for details of the spectral 
extractions). We also show the ‘host-only’ spectrum, which repre- 
sents an estimate of the host galaxy spectrum excluding the region 
of the GRB. Comparing the spectrum at the position of the GRB with 
the host-only spectrum clearly shows that certain lines, mostly from 
molecular H,, are much stronger at the position of GRB 221009A. 

We measured the global host properties by fittingthe global host 
spectrum, as well as our NIRCam photometry and HST photometry 
from ref. 31, using the stellar population modelling code Prospector^* 
(see Methods for details of the modelling procedure). The best- 
fit model spectrum and photometry compared with the observed 
data are shown in Fig. 6. We find that the host has a stellar mass of 
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log(M/Ms) = 9.61*027 and low stellar and gas-phase metallicities of 


log(Z,./Z) = —0.81*0 84 and log(Z;4,/Zo) = —0.96*005, respectively. 
This is one of the lowest metallicity environments of any LGRB, which 
is a class of objects that prefer low-metallicity galaxies? and it is, to 
our knowledge, the lowest metallicity environment of a GRB-SN to date. 
This may suggest that very low metallicity is required to produce a very 
energetic GRB. In addition, the galaxy exhibits a recent star formation 
rate (SFR) of SER, myr = 0.17 M, yr”. We also find that the galaxy exhibits 


a visual extinction of Ay = 0.6701 mag. This is consistent with our 


extinction constraints from the early-phaseJWST data (Methods) where 


we found a best-fit total extinction of Ay = 4.63+943 mag, which is in 


good agreement, within uncertainties, with the nominal Milky Way value 
plus the host galaxy extinction found here. Our SFR and host extinction 
values are consistent with those measured from Ha and Paa detected 
inanearly-phase X-shooter spectrum of GRB 221009A (ref. 20). 

We additionally model the spectrum at the site of the GRB, and 


find a similar gas-phase metallicity of log(Z;4,/Zo) = —0.94*011 and 


a lower stellar metallicity of log(Z./Za) = —1.66*026 compared with 


the global host galaxy (Extended Data Fig. 6), which indicates that 
the progenitor of GRB 2210094 originated from a low-metallicity 
environment. 


Strong H, emission at the explosion site 
We observe many narrow H, vibrational and rotational emission lines 
that appear strongest at the site of GRB 2210094, as highlighted in 
Fig. 5. Molecular hydrogen traces dense star-forming regions, con- 
sistent with a birth cloud of a massive stellar progenitor of a LGRB. 
Neglecting the afterglow itself, H, can be excited by both shocks (driven 
by, for example, stellar winds or Herbig-Haro objects) or directly by 
fluorescence”, Following ref. 61, we compared the ratios of H, lines 
in the range -1.1-2.1 um with various models of fluorescence versus 
collisionally excited emission from ref. 59 using a simple chi-squared 
metric (with appropriately propagated uncertainties). Owing to the 
strong detection of many lines (which are predicted to be absent inthe 
case of collisional excitation), we find a better match to fluorescence 
models, which is consistent with the dominant excitation method in 
many low-metallicity, blue compact dwarf galaxies”. Our measured line 
ratios and predicted model ratios are given in Extended Data Table 2. 
Only one other LGRB host, that of GRB 031203 (a relatively faint 
LGRB), has had a marginal detection of H, emission?. Molecular H in 
absorption due to vibrational excitation has also been observed in a 
small number of events (see, for example, refs. 63,64). Statistical studies 
of GRB hosts have found that most lack vibrationally excited H, (for 
example, ref. 65), which suggests that molecular H production is sup- 
pressed in LGRB hosts. It has been suggested that this suppression may 
be partially due to the low metallicities of the hosts” or the ongoing star 
formation, leading to a strong ionizing field”. The low metallicity and 
modest SFR measured by Prospector suggests that the latter may lead 
to observable H, emission in this event. These observations highlight 
the importance of the unique sensitivity and spatial resolution of the 
JWST when analysing the local environments of LGRB progenitors. 


Conclusions 

We present the detection with the JWST of an SN associated with the 
highly energetic event GRB 221009A. Despite being associated with 
the brightest GRB ever observed, the SN produced a modest amount 
(=0.09 M,) of radioactive Ni with no obvious signs of r-process nucleo- 
synthesis. The host galaxy suggests a very low-metallicity progeni- 
tor system—one of the lowest metallicity environments of all known 
LGRBs. In addition, the exceptional sensitivity and spatial resolution 
of the JWST allows us to detect a series of multiple molecular H, emis- 
sion lines at the position of the GRB, which is an observation that has 
been long anticipated. A secondary site of r-process nucleosynthesis 


remains an open question, which can observationally be uniquely 
probed by late-time IR spectroscopy. Our findings motivate future 
JWST campaigns to examine the nebular-phase spectra of supernovae 
associated with LGRBs. 


Methods 

Imaging observations and photometry 

We obtained imaging of GRB 221009A with the NIRCam using the 
F115W, F200W, F277W and F444W filters on 22 April 2023 starting 
at 07:08 UT. Each observation consisted of four dithered exposures 
with a total exposure time of 558 seconds. We downloaded the stage 3 
pipeline products from the Mikulski Archive for Space Telescopes 
(MAST) for analysis. GRB 2210094 is clearly detected, along with its 
host galaxy. 

To measure the flux from GRB 2210094 in each filter, we first mod- 
elled the host galaxy contribution using the galaxy profile fitting code 
GALFIT®®. We modelled the host galaxy as a single Sérsic component. 
During the fit, we masked the pixels containing the light from GRB 
2210094; we fitted for transient flux in a later step. The input, best-fit 
model and residual (best-fit model subtracted off) images are shown 
in Fig. 1. The residual image for the F115W filter shows no structure 
indicating that the galaxy light is well described by a single Sérsic 
component. The residual images in the three redder filters, however, 
exhibit remaining diffuse structure not captured by the model near the 
centre of the galaxy and to the northeast. Although there is no obvious 
evidence for such diffuse structure emanating from the position of GRB 
2210094, itis plausible that GRB 221009A is co-located with a brighter 
region of its host galaxy that is not captured by our galaxy model. Such 
a determination can only be made when the transient fades. 

Next, we performed point spread function (PSF) photometry on 
the residual images at the location of GRB 221009A. As WebbPSF only 
generates PSF models for use with stage 2 imaging data, we used the 
following custom procedure to generate drizzled PSFs for use with 
stage 3 data. We generated stage 2 images with model PSFs planted 
at the location of GRB 221009A and then ran these images through 
the stage 3 pipeline. We then used the drizzled PSF models for the 
PSF fitting of GRB 221009A. We found the following AB magnitudes 
in each filter: Mpusw = 25.10 + 0.05 mag, Mp ro9w = 24.12 + 0.11 mag, 
Meo77w = 23.77 + 0.05 mag and Mș444w = 23.22 + 0.08 mag (not corrected 
for extinction, which we assess in detail in Methods section ‘Constraints 
on foreground dust from early-time spectroscopy’). The uncertainties 
include the systematic uncertainty associated with the GALFIT model- 
ling procedure, which we estimated by comparing PSF photometry 
with and without galaxy subtraction with GALFIT. 


Spectroscopic observations and data reduction 
Late-time NIRSpec observations. We obtained spectra of GRB 
221009A on 20 April 2023 with the NIRSpec? onboard JWST (pro- 
gramme 2784; principal investigator (PI), Blanchard). Our observations 
began at 14:40 UT, corresponding to a rest-frame phase of 167.7 days 
since the Fermi GBM trigger. Spectra were taken with the S200A1 fixed 
slit and the medium-resolution gratings GI40M/F100LP and G235M/ 
F170LP, yielding wavelength coverage in the range -1-3 um. For each 
grating and filter set-up, we used five primary dithers and a total expo- 
sure time of 10,942 seconds. Owing to the small offset of GRB 221009A 
fromits host galaxy”, target acquisition was performed using an offset 
star to ensure proper centring of the source in the slit. 

We downloaded and inspected the pipeline products available on 
the MAST. A resolved trace is clearly present in the individual stage 2 
exposures and final combined stage 3 products, indicating a substan- 
tial host galaxy contribution. In addition, a compact trace spanning 
~2 pixels is apparent at the red end of the GI40M/F1OOLP spectrum 
and the G235M/F170LP spectrum at the expected location of GRB 
221009A within the slit. This trace is also at a consistent offset from 
thebrightest part ofthe resolved trace representing the centre ofthe 
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host galaxy, which confirms that this unresolved trace is the spectrum 
of GRB 221009A. 

The pipeline products available on MAST were reduced using 
nod-subtraction, the default background subtraction method for 
apointsource with multiple dithered exposures. Owing to the resolved 
nature of the overall trace, we re-reduced the data using the JWST 
Science Calibration Pipeline with nod-subtraction turned off to 
reduce the effect of subtracting source flux from itself. We then 
extract the one-dimensional spectrum of GRB 221009A from 
our re-reduced stage 3 combined and rectified two-dimensional (2D) 
spectra for further analysis. The final 2D spectra are shown in Extended 
Data Figs. 1 and 2 (see Fig. 1 for the slit orientation). 

We used the following extraction procedure to isolate the flux 
associated with GRB 221009A from the light of its host galaxy. We 
modelled the spatial profile of the overall trace as a two-component 
Gaussian with centres fixed at the position of GRB 221009A and the 
centre of the host galaxy. We fitted the total spatial profile, summed 
over all wavelengths, to determine the best-fit Gaussian width of each 
component. We then fitted this model, with widths fixed at these values, 
tothe spatial profile at each wavelength. We also fitted for a linear back- 
ground component determined from background regions located on 
both sides of the trace. The sum of the flux in each fitted Gaussian com- 
ponent thus represents the flux from GRB 221009A and its host galaxy 
as a function of wavelength. The resulting spectra of GRB 221009A in 
the G140M and G235M gratings are shown in Extended Data Figs. land 
2, respectively. We also show corresponding error spectra calculated 
from the 2D pipeline-generated error arrays for each grating. The fit 
to the host galaxy Gaussian component yields a ‘host-only’ spectrum. 

We note that the background exhibits evidence for PSF artefacts 
that are potentially from a nearby bright star (the pseudoperiodic 
signal at pixel rows -18-24 in the 2D frames; Extended Data Figs. 1 
and 2), which is the most likely explanation given the crowded nature 
ofthe field. Owing to the difficulty of accurately modelling this com- 
ponent of the background, our background regions exclude those 
containing such artefacts. Our fitted background therefore represents 
thesmooth underlying sky background. This may mean that the back- 
ground atthe location of GRB 221009A and its host galaxy is underesti- 
mated. However, we extract regions ofthe background containing the 
suspected PSF artefacts from a nearby star and find no evidence that 
these features are present in our extracted spectrum of GRB 221009A. 
In addition, the flux from these features decreases towards the spatial 
location of the GRB spectral trace. 


Combined G140M * G235M spectrum compared with photometry. 
As the photometry was obtained only two days after our NIRSpec 
spectra, we used the photometry to check the flux calibration of the 
spectra. We find that the fluxes in the four NIRCam filters are an excel- 
lent match to the flux calibration of the GI40M and G235M spectra. In 
addition, the two spectra agree in the wavelength region where the 
gratings overlap. In Fig. 2, we show the fluxes in each filter compared 
with the combined G140M + G235M spectrum. 


Host galaxy spectral extractions. As seen in Extended Data Figs. 1 
and 2, the host galaxy is resolved in our JWST/NIRSpec observations, 
extending across approximately ten rows in our 2D spectra and with 
numerous narrow emission lines. To study the global host properties, 
we extracted the entire trace including the position of the GRB. We 
note that there is significant variation of the strength of some emis- 
sion lines across the spatial extent of the galaxy, where several lines 
are stronger at the position of GRB 2210094. To identify these lines and 
assess any potential variation in the galaxy properties at the position 
of the GRB, we extracted a narrow aperture centred on the position of 
GRB 2210094. This differs from the Gaussian decomposition proce- 
dure described in Methods section ‘Late-time NIRSpec observations’ 
used to isolate the GRB spectrum, as here we are not modelling and 


subtracting theunderlying host spectrum; the goal hereisto measure 
the host properties at the position of the GRB. 


Archival NIRSpec/MIRI observations. We obtained archival spec- 
troscopic observations of GRB 221009A from JWST, observed with 
NIRSpec and the MIRI on 22 Oct 2022 (programme 2782; PI, Levan and 
originally presented in ref. 31). These observations correspond to 13.16 
and 13.2 days post burst, respectively. 

At this epoch, the NIRSpec observations were taken in the 
low-resolution PRISM mode, with spectral coverage from -0.5 jum 
to -5.5 um. The pipeline products from MAST reveal a clear, high 
signal-to-noise ratio trace in the 2D spectrum. The stage 3, reduced 
spectrum is consistent with that published in ref. 31, and we thus use 
itfor analysis inthis work without additional reductions. 

TheMIRIspectrum wastaken in the low-resolution spectroscopy 
mode with the P750L disperser. The automatic reduction of the MIRI 
spectrum failed, which was likely to be due to improper selection of the 
afterglow trace. We used the official MIRI reduction pipeline to manu- 
ally extract the spectrum from the stage 2 products, carefully selecting 
the correct trace and appropriate background from the nodded 2D 
image. At this epoch, the afterglow trace is clearly identified inthe 2D 
spectrum and easily isolated using a simple boxcar extraction. We note 
that MIRIis uncalibrated belowA x 4.5 um at the time of analysis, and we 
therefore removed data below this wavelength of the spectrum from 
analysis. The MIRI observations are qualitatively consistent with those 
of ref. 31 and are well matched to their near-simultaneous photometric 
observation in F560W. 


ALMA observations 

Following the seven epochs of ALMA observations of GRB 221009A 
through programme 2022.1.01433.T (PI, Laskar), we obtained two 
additional epochs with the same programme on 01 March 2023 at a 
mean time of 15:41 UT and on 11 April 2023 at a mean time of 07:55 UT, 
corresponding to 143.6514 and 183.7838 days in the observer frame, 
respectively (4124.80 and 7159.67 days in the rest frame). Both 
observations utilized two 4 GHz wide base-bands centred at 91.5 
and 103.5 GHz, respectively with J1924 — 2914 as bandpass calibra- 
tor and JI914 + 1636 as complex gain calibrator. The millimetre-band 
afterglow previously reported in ref. 25 was clearly detected in the 
pipeline-processed science-ready data products in the first of the 
two epochs and more weakly («4.70) detected in the second epoch. 
We performed photometry using imfit in the Common Astronomy 
Software Applications” and find a best-fit flux density in the two epochs 
of (163 + 22) wy and (104 + 23) wy (including a 5% systematic flux 
calibration uncertainty) at a mean frequency of 97.5 GHz, along with 
a position of RA - 19 h 13 m 03.50 s and dec.- 19°46’24.3” with an 
uncertainty of 0.1 in each coordinate (consistent across both epochs). 
Together with the last ALMA 97.5 GHz measurement reported in 
ref. 25, the temporal decline rate of the millimetre-band afterglow at 
-99-188 days after the burst (observer frame) is a, = 71.54 + 0.08, 
which implies an extrapolated millimetre-band flux density of 
(99 + 23) wy at the time of the NIRSpec observations (194 days, 
observer frame). 


Swift/XRT observations 

We downloaded the count-rate light curve of the X-ray afterglow 
of GRB 221009A from the Swift/XRT website”. Using the spectral 
parameters presented in ref. 25 (Milky Way (MW) absorption of 
Namw = 5.36 x 10? cm”, intrinsic absorption of Nyin = 1.35 x 10? cm? 
and photon index of [Fx = 1.8566), we converted the observed count rate 
to a flux density (F,) at 1 keV and obtained F, = (9.3 + 3.5) x 10? py 
at 196*55 days (observer frame; corresponding to 170.4*4? days, 
rest frame) after the burst. Comparing this with our ALMA obser- 
vations, we find that the spectral index between the ALMA and XRT 
observations at the time of the NIRSpec observations (2194 days, 
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observer frame) is Baima-xet = 0.63 + 0.03. We refer to this number 
elsewhere in the text as the ALMA-XRT power law, anchored to the 
inferred ALMA flux density of «0.1 mJy at the time of the NIRSpec 
observations. 


Constraints on foreground dust from early-time spectroscopy 
Given the location of GRB 2210094 in the Galactic plane (b = 4), we 
expected substantial extinction due to interstellar dust in the MW. 
Reference 36 estimates the MW extinction contribution to be 
Ay = 4.10 + 0.06, assuming the standard extinction factor Ry = 3.1. As 
noted in ref. 72, these dust maps can be unreliable for low galactic 
latitudes”. Furthermore, this measurement neglects host contribution; 
however, given the relatively low redshift, we do not expect to easily 
distinguish between dust arising from either the MW or host galaxy 
(see, for example, the results of ref. 20). For simplicity, we neglect 
redshift dependence of dust. 

Given the significant uncertainties expected, we opted to use 
the first epoch of NIRSpec/MIRI data to determine the appropriate 
extinction correction. We assumed that the spectrum is dominated by 
some unknown combination of an afterglow (power-law model) and 
athermal, SN-like component. Unless r-process material is mixed sig- 
nificantly within the ejecta, we do not expect ared thermal component 
at early times. As such, we assumed that the event is dominated by a 
power-law afterglow atA,,, 2 3 jum; below this wavelength, it is reason- 
able that a SN 1998bw-like event could contribute significant flux. We 
explore the systematic uncertainties associated with the extinction 
laws and assumptions on the SN contribution. 

Few extinction laws are calibrated across the full wavelength 
range covered by the NIRSpec/MIRI observations; a new extinction 
law describing A(A)/A(V) as a function of R(V) in the range -0.1-30 um 
has been recently presented”. We contrast this solution with the com- 
monly used extinction law described in ref. 37 to quantify systematic 
uncertainty from assumptions of dust laws. Given a prescribed dust 
law, we simultaneously fitted the observed day 13.2 (observer frame) 
spectrum to a power law (F, ~ v?) and extinction parameters A, and Ry 
using a Markov chain Monte Carlo sampler implemented in emcee”. 
Our models have four free parameters: the overall power-law normali- 
zation (‘norm’), the power-law index f, the dust A, and A, values, and 
a white noise scatter term. The scatter quantifies the uncertainty in 
JWST flux estimates as a fraction ofthe flux. We assumed a wide uniform 
prior for all parameters except normalization, in which we assumed a 
log-uniform prior. 

We first fitted using the dust law presented in ref. 38. Fitting all 
observed wavelengths A,,, < 8 um, we find £ = 0.39 + 0.01, A, = 4.37 + 
0.05 and Ry = 3.077095. At A < 2 um, we find that the residuals are 
consistent with O, suggesting no contribution from an additional ther- 
mal component. We next excluded wavelengths <2 um in the fitting 
process to test the possibility of contamination from either a SN-like 
or r-process thermal event. We find that when excluding these wave- 
lengths, the afterglow model overestimates the blue flux. 

Next, we fitted using the extinction law described in ref. 37 (that 
is, following the original analysis of ref. 31). Again, we emphasize that 
this extinction law is not calibrated for IR observations and simply 
extrapolates at these wavelengths. We simultaneously fitted the 
observed spectrum (A,, < 8 jum) to a power law and extinction model. 
We find £ = 0.41 + 0.01, Ay = 4.63*0 and Ry = 4.24*07^. This is signifi- 
cantly different (>30) from the results presented in ref. 31 when only 
accounting for statistical uncertainties, which we attribute to a tight 
prior (versus our flat prior) set by those authors. 

We report the results of our fits in Extended Data Table 1 and show 
these data, models and associated residuals in Extended Data Fig. 3. The 
residuals of both dust models show significant structure throughout 
the spectrum. We specifically compared the residuals to a spectrum of 
SN 1998bw taken 12 days post burst and scaled to the redshift of GRB 
2210094. We note that the statistical uncertainties and systematic 


difference between these two dust extinction models mean that we 
are unable to make a conclusive statement on the SN emission from 
theearly-timeJWST spectrum. This is a different conclusion from that 
of ref. 31, who, given their small statistical uncertainties, rule out SN 
1998bw-like thermal emission at early times. 


Constraints on the afterglow contribution 
Initial comparisons with previous supernovae. In Extended Data 
Fig. 5, we show our extinction-corrected spectrum (using the law in 
ref. 37 and best-fit parameters listed in Extended Data Table 1) com- 
pared with spectra of SN 1998bw^, the canonical SN Ic-BL associated 
with a GRB and SN 2013ge™, one of the few supernovae Ic with high S/N 
late-time NIR spectra, taken at +51 and +118 days after peak, respec- 
tively. To achieve complete overlap with the blue end of our spectrum, 
we combined the +51 day NIR spectrum of SN 1998bw with an optical 
spectrum taken at +73 days. We scaled the spectra of SN1998bw and SN 
2013ge to the distance of GRB 221009A and used their light curves*"? 
to normalize to their brightnesses at the phase of our GRB 221009A 
spectrum. Our spectrum of GRB 2210094 is brighter than the com- 
parison supernovae would be and relatively featureless with a different 
overall spectral shape, which is consistent with significant contami- 
nation from the afterglow. Our spectrum exhibits flux increasing at 
Az1.5 um, whereas the comparison supernovae exhibit declining flux. 
The emission features shown in Extended Data Fig. 4 exhibit similar, 
although slightly narrower, widths than the corresponding features in 
SN1998bw, SN 2013ge and SN 2014ad. Owing to the lack ofa late-time 
light curve for SN 2014ad, we scaled its spectrum to roughly match SN 
1998bw for comparison purposes. In addition, the lines in our JWST 
spectrum are diluted in strength and exhibit a different flux ratio. 
This, combined with the rising flux to the red, means that there is no 
simple luminosity scaling that will bring our spectrum of GRB 221009A 
into agreement with the comparison spectra. These observations are 
consistent with afterglow contamination. Furthermore, the lack of 
many strong SN features (for example, the strong P-Cygni features near 
=] um, 1.5 um and =2 um commonly seen in supernovae”) other than 
the two identified (Ca 11 NIR triplet and O rat «0.86 um and «0.92 pm, 
respectively) indicates that the SN associated with GRB 221009A is not 
substantially brighter than SN 1998bw and SN 2013ge. 


Constraints from contemporaneous ALMA and Swift/XRT observa- 
tions. Determining the afterglow contribution is critical for constrain- 
ing the presence of SN emission and a possible contribution from 
r-process material. First, we considered the power law formed by the 
ALMA and XRT observations that we obtained around the same phase as 
ourJWST observations. We analysed the residual spectrum by subtract- 
ing off the ALMA-XRT power law from our spectrum of GRB 221009A 
which we show, compared with SN 1998bw and SN 2013ge, in Extended 
Data Fig. 5. Although the resulting spectrum matches more closely the 
shape of the supernovae compared to the unsubtracted spectrum, in 
particular at the blue end, the shape at A z 1.5 pm still exhibits rising 
flux substantially different from the supernovae. Given the lack of 
strong emission features in this region, the most likely explanation is 
that the ALMA-XRT power-law model does not adequately capture the 
afterglow contribution. In section ‘No signs of r-process enrichment’ 
and Methods section ‘Comparison with r-process light curve models’, 
we consider whether this red excess could be due to emission from 
r-process material. 


Varying the afterglow contribution. Next we considered the best-fit 
power law from fitting our spectrum at Az 1.5 um (shown in Fig. 3) 
and analysed how the implied SN component changes with different 
afterglow normalizations. We scaled the best-fit power law by factors 
of 0.3, 0.6, 0.9 and 1.0 to generate four potential afterglow models, 
subtracted them from the spectrum and compared the resulting resi- 
dual spectra with SN 1998bw and SN 2013ge. In Supplementary Fig. 1, 
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we show the residual spectra and afterglow models for the four scalings. 
When scaled by 0.3 and 0.6, the residual spectra still exhibit flux rising 
tothe red, as inthe unsubtracted spectrum, which indicates that these 
models are not likely to account for all of the afterglow flux. 

In addition, there is a mismatch between the flux ratios of the 
expected emission lines. In other words, the detection of the Ca 11 NIR 
triplet at the strength we see, would imply the detection of other lines 
at strengths that are not observed. Of course, this reasoning relies on 
the assumption that the SN associated with GRB 221009A should appear 
similar to previous supernovae Ic/Ic-BL. Indeed it is possible that this 
SN may not show the same features as previous events and potentially 
an additional component from r-process emission which we assess in 
section ‘No signs of r-process enrichment’ and Methods section 'Com- 
parison with r-process light curve models’. However, the lack of strong 
lines in this region indicates that the SN associated with GRB 221009A 
is likely to be fainter than these models suggest and the afterglow is 
correspondingly brighter (as found when performing a joint SN + after- 
glow fit; see section ‘Isolating the SN signal’), such that most emission 
lines are diluted with respect to the continuum andare not detectable. 

Ifinstead the best-fit power law is scaled by 0.9, the residual spec- 
trum appears consistent with the comparison spectra and is a close 
match to the overall flux level of SN 1998bw. Note that this is similar to 
the best-fit scaling (0.94) when performing the joint SN + afterglow 
fit as described in ‘Isolating the SN signal’ and shown in Fig. 3. Larger 
afterglow contributions (for example, scaling by 1.0) yield an overall 
steeper slope, inconsistent with the comparison objects. 


Comparison with r-process light curve models 

We also considered the r-process enriched SN light curve models in 
ref. 52. In Supplementary Fig. 2, we show the J-H and J-K colour evolu- 
tion of these models, for an SN Ic-BL with a typical simulated ejecta 
mass of 3.96 M,, a °°Ni mass of 0.33 M,, an r-process material mass 
of 0.03 M, and various levels of r-process mixing, from no mixing to 
nearly fully mixed, compared with the colours of the SN component 
of GRB 2210094 under different afterglow assumptions. We calcu- 
latedJ-H andJ-K colours by convolving the filter bandpasses with our 
NIRSpec spectrum after subtracting the afterglow models. We show 
the resulting colours for the afterglow models considered in section 
‘Constraints on the afterglow contribution’ (the ALMA-XRT power law 
and the best-fit power law from fitting the red end of our spectrum with 
various normalizations; Extended Data Fig. 5 and Supplementary Fig. 1). 

The J-K colours of the afterglow-subtracted spectra match the 
r-process enriched models when scaling the best-fit power law by 
<0.9. Decreasing the afterglow contribution leads to more residual red 
light, leading to redder colours. When scaling by <0.6, including the 
ALMA-XRT model, the J-K colours, if reddened due to r-process mate- 
rial, would imply significant mixing. In this case, strong broad emission 
lines from r-process elements would be expected, as seen in the MHD 
model in Fig. 4 but not in our data. In addition, for a given afterglow 
contribution, the J-H colours imply a different degree of r-process 
mixing than the J-K colours, which suggests that the reddening source 
is not dueto r-process emission. 

In Supplementary Fig. 2, we also show the colours of SN 1998bw 
and SN2013ge calculated fromtheir late-time NIR spectra. SN 1998bw 
is notably blue—bluer even than the models without r-process—which 
suggests that these models do not fully capture the range of possible 
spectral energy distributions oftypical GRB supernovae. SN 2013ge is 
notably red, whichis consistent with the r-process enriched models for 
a mixing fraction of ~10%. This event, however, exhibits a clear exam- 
ple of carbon monoxide emission increasing the flux in the K-band. 
These comparisons highlight that, without spectra, other sources 
of reddening are difficult to disentangle from that due to r-process 
material. Similar conclusions have been drawn from studies of large 
samples of supernovae Ic-BL light curves?. We note that the spectrum 
of GRB 221009A after subtracting the best-fit power law scaled by 0.9, 


which yields a good visual match to SN 1998bw and SN 2013ge (Sup- 
plementary Fig. 1), exhibits a J-K colour that is -0.1 mag redder than 
the no r-process model. However, as can be seen in Supplementary 
Fig. 1, there is an upturn in the spectrum in the K-band at the expected 
location of the first overtone carbon monoxide emission, similar to 
that seen in SN 2013ge. 


Host galaxy modelling 

We used Prospector™, a Bayesian galaxy spectral energy distribution 
(SED) fitting code to simultaneously fit the global host galaxy photom- 
etry and spectroscopy. Additionally, we fitted the spectrum extracted 
at the position of the GRB to compare the global host properties and 
those at the GRB position. We adopted the MIST isochrones” and the 
C3K stellar spectral libraries in the Flexible Stellar Population Syn- 
thesis?" framework. The stellar population is described by redshift, 
stellar mass, velocity dispersion, stellar metallicity and a step function 
non-parametric star formation history with 14 time bins”. The nebular 
emission was parametrized by gas-phase metallicity and ionization 
parameter using the CLOUDY grid in ref. 79. We simultaneously fit- 
ted simple Gaussians to lines that are not included in our emission 
line model that assumes that all the emission is powered by the stars, 
namely, the He I, [Fe r1] and H, emission lines, with the same kinemat- 
ics but free amplitudes as our CLOUDY grid. We assumed a flexible 
two-component dust attenuation model that accounted for birth cloud 
and diffuse dust separately??. Variation in the shape of the attenuation 
curve was enabled using a power-law modification to a Calzetti curve?'. 
Wealso incorporated the contribution of dust emission to the infrared 
photometry using a three-parameter model? To fit the spectros- 
copy and photometry together, we marginalized over the shape of the 
observed spectrum (thereby avoiding any wavelength-dependent flux 
calibration issues) with a polynomial; in this manner, the normalization 
and shape of the SED was entirely determined by the photometry, or not 
constrained at all for the GRB position, where there is no photometry. 
We also included a jitter parameter that inflated the spectroscopy 
uncertainties to account for imperfect JWST flux calibration and slit 
losses, and found typical values of 1.5-2, which are consistent with 
other early JWST spectroscopic analyses®*. Finally, we used a pixel 
outlier model to downweight pixels that were not consistent with our 
model®, which were typically identified at a 1-296 level. In summary, the 
SED modelfor the host galaxy fit has 28 free parameters, and the fit to 
thespectrum at the GRB position has 24 free parameters. 


Data availability 
TheJWST data analysed in this work associated with programmes 2784 
and 2782 are publicly available on the MAST archive. 


Code availability 

The software tools used in this work (JWST Science Calibration Pipe- 
line, GALFIT, WebbPSF, Prospector and Common Astronomy Software 
Applications) are publicly available. 


References 

1. Cowan, J. J. et al. Origin of the heaviest elements: the rapid 
neutron-capture process. Rev. Mod. Phys. 93, 015002 (2021). 

2. Siegel, D. M. r-process nucleosynthesis in gravitational-wave and 
other explosive astrophysical events. Nat. Rev. Phys. 4, 306-318 
(2022). 

3. Symbalisty, E. & Schramm, D. N. Neutron star collisions and the 
r-process. Astrophys. Lett. 22, 143-145 (1982). 

4. Eichler, D., Livio, M., Piran, T. & Schramm, D. N. Nucleosynthesis, 
neutrino bursts and y-rays from coalescing neutron stars. Nature 
340, 126-128 (1989). 

5. Arcavi, l. et al. Optical emission from a kilonova following a 
gravitational-wave-detected neutron-star merger. Nature 551, 
64-66 (2017). 


Nature Astronomy 


Article 


https://doi.org/10.1038/s41550-024-02237-4 


10. 


12. 


13. 


14. 


15. 


16. 


17. 


18. 


19. 


20. 


21. 


22. 


23. 


24. 


25. 


26. 


27. 


28. 


29. 


Cowperthwaite, P. S. et al. The electromagnetic counterpart 

of the binary neutron star merger LIGO/Virgo GW170817. II. UV, 
optical, and near-infrared light curves and comparison to kilonova 
models. Astrophys. J. Lett. 848, L17 (2017). 

Chornock, R. et al. The electromagnetic counterpart of the 
binary neutron star merger LIGO/Virgo GW170817. IV. Detection 
of near-infrared signatures of r-process nucleosynthesis with 
Gemini-South. Astrophys. J. Lett. 848, L19 (2017). 

Drout, M. R. et al. Light curves of the neutron star merger 
9w170817/sss17a: implications for r-process nucleosynthesis. 
Science 358, 1570-1574 (2017). 

Smartt, S. et al. A kilonova as the electromagnetic counterpart to 
a gravitational-wave source. Nature 551, 75-79 (2017). 

Côté, B. et al. Advanced LIGO constraints on neutron star mergers 
and r-process sites. Astrophys. J. 836, 230 (2017). 

Naiman, J. P. et al. First results from the illustrisTNG simulations: 
atale of two elements - chemical evolution of magnesium and 
europium. Mon. Not. R. Astron. Soc. 477, 1206-1224 (2018). 

Cóté, B. et al. Neutron star mergers might not be the only source 
of r-process elements in the Milky Way. Astrophys. J. 875, 106 
(2019). 

Hotokezaka, K., Beniamini, P. & Piran, T. Neutron star mergers as 
sites of r-process nucleosynthesis and short gamma-ray bursts. 
Int. J. Mod. Phys. D 27, 1842005 (2018). 

Kobayashi, C., Karakas, A. |. & Lugaro, M. The origin of elements 
from carbon to uranium. Astrophys. J. 900, 179 (2020). 

Siegel, D. M., Barnes, J. & Metzger, B. D. Collapsars as a major 
source of r-process elements. Nature 569, 241-244 (2019). 

Just, O., Aloy, M., Obergaulinger, M. & Nagataki, S. r-process 
viable outflows are suppressed in global alpha-viscosity models 
of collapsar disks. Astrophys. J. Lett. 934, L30 (2022). 

Burns, E. et al. Grb 221009a: the boat. Astrophys. J. Lett. 946, L31 
(2023). 

Lesage, S. et al. Fermi-GBM discovery of GRB 2210094: an 
extraordinarily bright GRB from onset to afterglow. Astrophys. J. 
Lett. 952, LA2 (2023). 

Williams, M. A. et al. GRB 2210094: discovery of an exceptionally 
rare nearby and energetic gamma-ray burst. Astrophys. J. Lett. 
946, L24 (2023). 

Malesani, D. B. et al. The brightest GRB ever detected: GRB 
221009A as a highly luminous event at z=0.151. Preprint at 
arXiv:2302.07891 (2023). 

Woosley, S. E. Gamma-ray bursts from stellar mass accretion disks 
around black holes. Astrophys. J. 405, 273-277 (1993). 

Fryer, C. L., Woosley, S. E. & Hartmann, D. H. Formation rates of 
black hole accretion disk gamma-ray bursts. Astrophys. J. 526, 
152-177 (1999). 

MacFadyen, A. |. & Woosley, S. E. Collapsars: gamma-ray bursts 
and explosions in ‘failed supernovae’. Astrophys. J. 524, 262-289 
(1999). 

Cano, Z., Wang, S.-Q., Dai, Z.-G. & Wu, X.-F. The observer's guide 
to the gamma-ray burst supernova connection. Adv. Astron. 
https://doi.org/10.1155/2017/8929054 (2017). 

Laskar, T. et al. The radio to GeV afterglow of GRB 221009A. 
Astrophys. J. Lett. 946, L23 (2023). 

O'Connor, B. et al. A structured jet explains the extreme GRB 
221009A. Sci. Adv. 9, eadi1405 (2023). 

Bright, J. S. et al. Precise measurements of self-absorbed rising 
reverse shock emission from gamma-ray burst 221009A. 

Nat. Astron. 7, 986-995 (2023). 

Shrestha, M. et al. Limit on supernova emission in the brightest 
gamma-ray burst, GRB 2210094. Astrophys. J. Lett. 946, L25 (2023). 
Fulton, M. D. et al. The optical light curve of GRB 2210094: the 
afterglow and the emerging supernova. Astrophys. J. Lett. 946, 
L22 (2023). 


30. 


31. 


32. 


33. 


34. 


35. 


36. 


37. 


38. 


39. 


40. 


41. 


42. 


43. 


44. 


45. 


46. 


47. 


48. 


49. 


50. 


51. 


52. 


53. 


54. 


Srinivasaragavan, G. P. et al. A sensitive search for supernova 
emission associated with the extremely energetic and nearby 
GRB 221009A. Astrophys. J. Lett. 949, L39 (2023). 

Levan, A. J. et al. The first JWST spectrum of a GRB afterglow: no 
bright supernova in observations of the brightest GRB of all time, 
GRB 221009A. Astrophys. J. Lett. 946, L28 (2023). 

Rastinejad, J. C. et al. A kilonova following a long-duration 
gamma-ray burst at 350 Mpc. Nature 612, 223-227 (2022). 
Levan, A. et al. Heavy-element production in a compact object 
merger observed by JWST. Nature 626, 737-741 (2024). 

Troja, E. et al. A nearby long gamma-ray burst from a merger of 
compact objects. Nature 612, 228-231 (2022). 

Yang, Y.-H. et al. A lanthanide-rich kilonova in the aftermath of a 
long gamma-ray burst. Nature 626, 742-745 (2024). 

Schlafly, E. F. & Finkbeiner, D. P. Measuring reddening with 
Sloan Digital Sky Survey stellar spectra and recalibrating SFD. 
Astrophys. J. 737, 103 (2011). 

Fitzpatrick, E. L. Correcting for the effects of interstellar 
extinction. Publ. Astron. Soc. Pac. 111, 63-75 (1999). 

Gordon, K. D. et al. One relation for all wavelengths: the far- 
ultraviolet to mid-infrared Milky Way spectroscopic R(V)- 
dependent dust extinction relationship. Astrophys. J. 950, 86 
(2023). 

Shahbandeh, M. et al. Carnegie Supernova Project-ll: 
near-infrared spectroscopy of stripped-envelope core-collapse 
supernovae. Astrophys. J. 925, 175 (2022). 

Patat, F. et al. The metamorphosis of SN 1998bw. Astrophys. J. 
555, 900-917 (2001). 

Clocchiatti, A., Suntzeff, N. B., Covarrubias, R. & Candia, P. The 
ultimate light curve of SN 1998bw/GRB 980425. Astron. J. 141, 163 
(2011). 

Drout, M. R. et al. The double-peaked SN 2013ge: a type Ib/c SN 
with an asymmetric mass ejection or an extended progenitor 
envelope. Astrophys. J. 821, 57 (2016). 

Iwamoto, K. et al. Ahypernova model for the supernova 
associated with the y-ray burst of 25 April 1998. Nature 395, 
672-674 (1998). 

Sollerman, J. et al. SN 1998bw at late phases. Astrophys. J. Lett. 
537, L127-L130 (2000). 

Nakamura, T., Mazzali, P. A., Nomoto, K. & lwamoto, K. Light curve 
and spectral models for the hypernova SN 1998BW associated 
with GRB 980425. Astrophys. J. 550, 991-999 (2001). 

Mazzali, P. A., Nomoto, K., Patat, F. & Maeda, K. The nebular 
spectra of the hypernova SN 1998bw and evidence for 
asymmetry. Astrophys. J. 559, 1047-1053 (2001). 

Maeda, K. et al. A two-component model for the light curves of 
hypernovae. Astrophys. J. 593, 931-940 (2003). 

Hjorth, J. The supernova-gamma-ray burst-jet connection. 
Philos. Trans. R. Soc. A 371, 20120275-20120275 (2013). 

Anand, S. et al. Collapsars as sites of r-process nucleosynthesis: 
systematic near-infrared follow-up of type Ic-BL supernovae. 
Astrophys. J. 962, 68 (2024). 

Siegel, D. M. et al. ‘Super-kilonovae’ from massive collapsars as 
signatures of black hole birth in the pair-instability mass gap. 
Astrophys. J. 941, 100 (2022). 

Barnes, J. & Duffell, P. C. Hydrodynamic mixing of accretion disk 
outflows in collapsars: implications for r-process signatures. 
Astrophys. J. 952, 96 (2023). 

Barnes, J. & Metzger, B. D. Signatures of r-process enrichment in 
supernovae from collapsars. Astrophys. J. Lett. 939, L29 (2022). 
Planck Collaboration. et al. Planck 2018 results. VI. Cosmological 
parameters. Astron. Astrophys. 641, A6 (2020). 

Johnson, B. D., Leja, J., Conroy, C. & Speagle, J. S. Stellar 
population inference with Prospector. Astrophys. J. Suppl. Ser. 
254, 22 (2021). 


Nature Astronomy 


Article 


https://doi.org/10.1038/s41550-024-02237-4 


55. Levesque, E. M., Kewley, L. J., Berger, E. & Zahid, H. J. The host 
galaxies of gamma-ray bursts. II. A mass-metallicity relation for 
long-duration gamma-ray burst host galaxies. Astron. J. 140, 
1557-1566 (2010). 

56. Graham, J. F. & Fruchter, A. S. The metal aversion of long-duration 
gamma-ray bursts. Astrophys. J. 774, 119 (2013). 

57. Perley, D. A., Niino, Y., Tanvir, N. R., Vergani, S. D. & Fynbo, J. P. U. 
Long-duration gamma-ray burst host galaxies in emission and 
absorption. Space Sci. Rev. 202, 111-142 (2016). 

58. Graham, J. F., Schady, P. & Fruchter, A. S. A surprising lack of 
metallicity evolution with redshift in the long gamma-ray burst 
host galaxy population. Preprint at arXiv:1904.02673 (2019). 

59. Black, J. H. & Van Dishoeck, E. F. Fluorescent excitation of 
interstellar H2. Astrophys. J. 322, 412-449 (1987). 

60. Vanzi, L. & Rieke, G. Infrared spectroscopy of blue dwarf galaxies. 
Astrophys. J. 479, 694-701 (1997). 

61. Izotov, Y. |. & Thuan, T. X. Near-infrared spectroscopy of five blue 
compact dwarf galaxies: Il Zw 40, Mrk 71, Mrk 930, Mrk 996, and 
SBS 0335- 052e. Astrophys. J. 734, 82 (2011). 

62. Wiersema, K. et al. Infrared molecular hydrogen lines in GRB host 
galaxies. Mon. Not. R. Astron. Soc. 481, 1126-1132 (2018). 

63. Prochaska, J. et al. The first positive detection of molecular gas in 
a GRB host galaxy. Astrophys. J. Lett. 691, L27-L32 (2009). 

64. Heintz, K. et al. New constraints on the physical conditions 
in H;-bearing GRB-host damped Lyman-a absorbers. Astron. 
Astrophys. 629, A131 (2019). 

65. Tumlinson, J., Prochaska, J. X., Chen, H.-W., Dessauges-Zavadsky, M. & 
Bloom, J. S. Missing molecular hydrogen and the physical 
conditions of GRB host galaxies. Astrophys. J. 668, 667-673 
(2007). 

66. Petitjean, P., Ledoux, C., Noterdaeme, P. & Srianand, R. Metallicity 
asacriterion to select H,-bearing damped Lyman-systems. 
Astron. Astrophys. 456, L9-L12 (2006). 

67. Whalen, D., Prochaska, J. X., Heger, A. & Tumlinson, J. The 
molecular hydrogen deficit in gamma-ray burst afterglows. 
Astrophys. J. 682, 1114-1123 (2008). 

68. Peng, C. Y., Ho, L. C., Impey, C. D. & Rix, H.-W. Detailed decompo- 
sition of galaxy images. II. Beyond axisymmetric models. Astron. J. 
139, 2097-2129 (2010). 

69. Jakobsen, P. et al. The Near-Infrared Spectrograph (NIRSpec) on 
the James Webb Space Telescope l. Overview of the instrument 
and its capabilities. Astron. Astrophys. 661, A80 (2022). 

70. CASA Team. et al. CASA, the Common Astronomy Software 
Applications for radio astronomy. Publ. Astron. Soc. Pac. 134, 
114501 (2022). 

7A. Swift/XRT light curves of GRB 221009A (UKSSDC, accessed 18 May 
2023); https://www.swift.ac.uk/xrt curves/01126853 

72. Kann, D. et al. GRANDMA and HXMT observations of GRB 
221009a: the standard-luminosity afterglow of a hyper-luminous 
gamma-ray burst. Astrophys. J. Lett. 948, L12 (2023). 

73. Popowski, P., Cook, K. H. & Becker, A. C. The large-scale extinction 
map of the galactic bulge from the MACHO project photometry. 
Astron. J. 126, 2910-2921 (2003). 

74. Foreman-Mackey, D., Hogg, D. W., Lang, D. & Goodman, J. mcee: 
the MCMC hammer. Publ. Astron. Soc. Pac. 125, 306-312 (2013). 

75. Choi, J. et al. Mesa Isochrones and Stellar Tracks (MIST). |. Solar- 
scaled models. Astrophys. J. 823, 102 (2016). 

76. Conroy, C., Gunn, J. E. & White, M. The propagation of uncertainties 
in stellar population synthesis modeling. |. The relevance of 
uncertain aspects of stellar evolution and the initial mass function 
to the derived physical properties of galaxies. Astrophys. J. 699, 
486-506 (2009). 

77. Conroy, C. & Gunn, J. E. The propagation of uncertainties in 
stellar population synthesis modeling. IIl. Model calibration, 
comparison, and evaluation. Astrophys. J. 712, 833-857 (2010). 


78. Leja, J., Carnall, A. C., Johnson, B. D., Conroy, C. & Speagle, J. S. 
How to measure galaxy star formation histories. Il. Nonparametric 
models. Astrophys. J. 876, 3 (2019). 

79. Byler, N., Dalcanton, J. J., Conroy, C. & Johnson, B. D. Nebular 
continuum and line emission in stellar population synthesis 
models. Astrophys. J. 840, 44 (2017). 

80. Kriek, M. & Conroy, C. The dust attenuation law in distant galaxies: 
evidence for variation with spectral type. Astrophys. J. Lett. 775, 
L16 (2013). 

81. Noll, S. et al. Analysis of galaxy spectral energy distributions 
from far-UV to far-IR with CIGALE: studying a SINGS test sample. 
Astron. Astrophys. 507, 1793-1813 (2009). 

82. Draine, B. T. & Li, A. Infrared emission from interstellar dust. IV. The 
silicate-graphite-PAH model in the post-Spitzer era. Astrophys. J. 
657, 810-837 (2007). 

83. Belli, S. et al. Massive and multiphase gas outflow in a quenching 
galaxy at z=2.445. Preprint at arXiv:2308.05795 (2023). 

84. Brinchmann, J. High-z galaxies with JWST and local analogues - 
it is not only star formation. Mon. Not. R. Astron. Soc. 525, 
2087-2106 (2023). 

85. de Graaff, A. et al. lonised gas kinematics and dynamical 
masses of z26 galaxies from JADES/NIRSpec high-resolution 
spectroscopy. Preprint at arXiv:2308.09742 (2023). 

86. Hogg, D. W., Bovy, J. & Lang, D. Data analysis recipes: fitting a 
model to data. Preprint at arXiv:1008.4686 (2010). 


Acknowledgements 

P.K.B. acknowledges support from a CIERA Postdoctoral Fellowship. 
V.AV. acknowledges support by the NSF through grant AST-2108676. 
The authors thank O. Fox and M. Shahbandeh for assistance in 
reducing the MIRI spectrum. This study was enabled in part by 

a Radboud Excellence Fellowship from Radboud University in 
Nijmegen, Netherlands. B.D.M. acknowledges support from the 
National Science Foundation (grant number AST-2002577). This work 
is based on observations made with the NASA/ESA/CSA James Webb 
Space Telescope. The data were obtained from the Mikulski Archive 
for Space Telescopes at the Space Telescope Science Institute, 

which is operated by the Association of Universities for Research 

in Astronomy, Inc., under NASA contract NAS 5-03127 for JWST. 

These observations are associated with programmes 2784 and 2782. 
Support for programme 2784 was provided by NASA through a grant 
awarded to P.K.B. from the Space Telescope Science Institute, which is 
operated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS 5-03127. This paper makes use of 

the following ALMA data: ADS/JAO.ALMA#2022.1.01433.T. ALMA is a 
partnership of ESO (representing its member states), NSF (USA) and 
NINS (Japan), together with NRC (Canada), MOST and ASIAA (Taiwan) 
and KASI (Republic of Korea), in cooperation with the Republic of 
Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and 
NAOJ. This work makes use of data supplied by the UK Swift Science 
Data Centre at the University of Leicester and of data obtained through 
the High Energy Astrophysics Science Archive Research Center 
On-line Service, provided by the NASA/Goddard Space Flight Center. 


Author contributions 

P.K.B. led the overall project from data proposal and acquisition to 
data analysis and the writing of the paper. P.K.B. is the PI of JWST 
programme 2784, the primary data analysed in this work. V.AV., as 
co-PI of JWST programme 2784, co-conceived the original JWST 
proposal. V.AV. contributed significant analysis of both our late-time 
JWST data and the archival JWST data from programme 2782 and 

a significant contribution to the writing of the paper. R.C., co-PI 

of JWST programme 2784, made significant contributions to the 
original JWST proposal and provided significant comments on the 
analysis and paper. T.L. analysed the ALMA and XRT data, wrote the 


Nature Astronomy 


Article 


https://doi.org/10.1038/s41550-024-02237-4 


Methods sections describing that data and provided comments on 
the paper. Y.L. and J.L. carried out the Prospector modelling of the 
host galaxy spectra and photometry and wrote the Methods section 
about such modelling. J.P. performed the PSF fitting of GRB 221009A 
in the NIRCam images. E.B. and R.M. provided substantial comments 
on the paper and the original JWST proposal. K.A., Y.C. and T.E. 
assisted with the acquisition of the ALMA and XRT data and provided 
comments on the paper. J.M.P. assisted with the reduction of the 
NIRSpec spectroscopy. D.S. provided the r-process spectral models 
for comparison and substantial comments on the paper. B.M., J.B., 
D.K., H.S., N.L. and S.K Y. provided comments on the paper and the 
original JWST proposal. A.R. assisted with the NIRCam photometry 
and commented on the paper. 


Competing interests 
The authors declare no competing interests. 


Additional information 
Extended data is available for this paper at 
https://doi.org/10.1038/s41550-024-02237-4. 


Supplementary information The online version contains supplementary 
material available at https://doi.org/10.1038/s41550-024-02237-4. 


Correspondence and requests for materials should be addressed to 
Peter K. Blanchard. 


Peer review information Nature Astronomy thanks Lauren Rhodes and 
the other, anonymous, reviewer(s) for their contribution to the peer 
review of this work. 


Reprints and permissions information is available at 
www.nature.com/reprints. 


Publisher's note Springer Nature remains neutral with regard 
to jurisdictional claims in published maps and institutional 
affiliations. 


Open Access This article is licensed under a Creative Commons 
Attribution 4.0 International License, which permits use, sharing, 
adaptation, distribution and reproduction in any medium or format, 
as long as you give appropriate credit to the original author(s) and the 
source, provide a link to the Creative Commons licence, and indicate 
if changes were made. The images or other third party material in this 
article are included in the article's Creative Commons licence, unless 
indicated otherwise in a credit line to the material. If material is not 
included in the article's Creative Commons licence and your intended 
use is not permitted by statutory regulation or exceeds the permitted 
use, you will need to obtain permission directly from the copyright 
holder. To view a copy of this licence, visit http://creativecommons. 
org/licenses/by/A.O/. 


€ The Author(s) 2024 


‘Center for Interdisciplinary Exploration and Research in Astrophysics (CIERA), Northwestern University, Evanston, IL, USA. ?Center for Astrophysics 

| Harvard & Smithsonian, Cambridge, MA, USA. ?Department of Astronomy, University of California, Berkeley, CA, USA. “Department of Physics & 
Astronomy, University of Utah, Salt Lake City, UT, USA. "Department of Astrophysics/IMAPP, Radboud University, Nijmegen, The Netherlands. Department 
of Astronomy & Astrophysics, The Pennsylvania State University, University Park, PA, USA. "Institute for Gravitation and the Cosmos, The Pennsylvania 
State University, University Park, PA, USA. ?Institute for Computational & Data Sciences, The Pennsylvania State University, University Park, PA, USA. 
?Space Telescope Science Institute, Baltimore, MD, USA. "Department of Physics, University of California, Berkeley, CA, USA. "Department of Astronomy/ 
Steward Observatory, Tucson, AZ, USA. "Kavli Institute for Theoretical Physics, University of California, Santa Barbara, CA, USA. “Department of Physics 
and Columbia Astrophysics Laboratory, Columbia University, New York, NY, USA. "Center for Computational Astrophysics, Flatiron Institute, New York, 
NY, USA. "Department of Physics and Astronomy, Northwestern University, Evanston, IL, USA. “Institute of Physics, University of Greifswald, Greifswald, 


Germany. "Department of Physics, University of Guelph, Guelph, Ontario, Canada. 


e-mail: peter.blanchard@northwestern.edu 


Nature Astronomy 


Article https://doi.org/10.1038/s41550-024-02237-4 


1.0 


0.8 


MEE trm npe decre] 


Pixel row 


1.35 1.4 1.45 1.5 1.55 
Wavelength (um) 


a rmm GRB221009A 
0.81 —— Error Spectrum 
ES 
3 0.6 
x 
30.4 
(rs 
0.2 
0.0 : . 
1.0 1.05 1.1 1.15 1.2 1.25 1:3 1.35 1.4 1.45 1:5 1.55 1.6 1.65 1:7 1.75 1.8 
Wavelength (um) 
Extended Data Fig. 1| NIRSpec/G140M Spectrum. Top: Final combined 2D is likely due to the diffraction spike of a nearby star. Spatially resolved emission 
NIRSpec/G140M spectrum of GRB 2210094 resulting from our re-reduction. The lines from the host galaxy are also detected. Bottom: ID spectrum of GRB 


trace is clearly dominated by the spatially resolved host galaxy. A broad emission 221009A extracted via a two-component Gaussian fit to the spatially resolved 
feature is visible near ~ lum (white arrow) at the expected spatial location of GRB trace to isolate the spectrum of GRB 2210094 from the host galaxy contribution. 


221009A (red arrow). The background below the trace contains structure that 
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host galaxy contribution. Notably, several host galaxy emission lines appear to be 
stronger at the position of GRB 221009A than the rest of the galaxy, resulting in 
significant excess flux from these lines appearing in the GRB 221009A spectrum 


Extended Data Fig. 2 | NIRSpec/G235M Spectrum. Top: Final combined 2D 
NIRSpec/G235M spectrum of GRB 221009A resulting from our re-reduction. 
The trace from GRB 221009A (red arrow) is more clearly visible at the red end. 
Bottom: 1D spectrum of GRB 2210094 extracted via a two-component Gaussian (see Fig. 5 for line identifications). 
fit to the spatially resolved trace to isolate the spectrum of GRB 221009A from the 
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Extended Data Fig. 3 | Extinction Fits to the Early JWST Data. Top: Comparison here, we fit A < 8m to avoid this feature. Bottom: Model residuals compared to 
of the early-time NIRSpec/PRISM and MIRI spectra (black) along with two models SN1998bw (black; [40]). Again, individual lines represent independent draws 


Residual F, [Jy] 


for the afterglow and dust law ([37] in blue, and [38] in orange). Each line of the from the model posterior. For both laws, strong systematic residuals are found 
model represents a draw from the posterior. Note that there is a silicate feature across the full wavelength range. An SN 1998bw-like supernova cannot be ruled 
at ~ 10um which is not properly modeled without a detailed dust composition; out. 
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Extended Data Fig. 4 | Broad SN Features. Left: Zoom-in on the blue end of our 
spectrum of GRB 2210094 highlighting the broad emission features we attribute 
to the Call NIR triplet and OI. Also shown are comparison spectra of SN 1998bw 
(orange) and SN 2013ge (gold) both scaled to the distance of GRB 221009A and 


their brightness at the phase of our observations, as well as SN 2014ad (magenta) 
arbitrarily scaled. Right: Spectrum of GRB 221009A after subtracting our best-fit 


afterglow model from the joint SN+afterglow fit described by the power law 
F, œ y 0762007. 
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Extended Data Fig. 5| Comparisons with Previous SNe and Assessment of the 
ALMA-XRT Power Law. Left: Our NIRSpec spectrum of GRB 2210094, corrected 
for extinction using the law from ref. 37 and best-fit parameters from the fitting 
described in the Methods (blue), compared to ground-based late-time NIR 
spectra of SN1998bw (orange) and SN 2013ge (gold) scaled to the distance of GRB 
221009A and the same phase of our observations. Right: Comparison between 
SN 1998bw and SN 2013ge and our spectrum of GRB 221009A (smoothed) after 


rw 
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subtracting an estimate of the afterglow contribution as described by the power 
lawconnecting ALMA and Swift/XRT observations taken around the same phase 
as our JWST data (F, « v ?5': black line). Significant flux rising toward the red 
remains in the resulting subtracted spectrum, inconsistent with the comparison 
SNe, indicating this model likely does not adequately describe the afterglow 
contribution at these wavelengths. 
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Extended DataFig. 6 | Galaxy Model Posteriors. Comparison ofthe Prospector model parameter posteriors for the galaxy spectra extracted at the location of GRB 


221009A (orange) and the total galaxy (blue). The gas-phase metallicities are similar, butthe GRB position exhibits a lower stellar metallicity. 
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Extended Data Table 1| Extinction Fitting Results 


Dust Law _log;,)(Norm, Jy) B Av Ry 
0.01 0.01 0.05 0.04 
[38] 3857991 0.397 9j 137 ros 3.07 gs 
[37] —3.83" 901 0.41759; 463766, 4.2476 64 


Best-fit extinction parameters and 1c uncertainties from fitting the +13.2 day (observer frame) NIRSpec/MIRI data. 
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Extended Data Table 2 | H, Emission Line Measurements 


Line (um) Flux / 107? (Jy) Measured Ratio Fluorescence Model Shock Model 
H» 1.233 2.15 «E 2.95 0.17 XE 0.24 0.47 0.01 
Hə 1.314 5.32 + 1.19 0.43 + 0.18 0.53 0.01 
Hə 1.957 20.70 + 1.43 1.66 + 0.62 NA NA 
Hə 2.033 11.16 + 0.92 0.89 + 0.33 0.56 0.37 
Hə 2.073 3.37 + 0.58 0.27 2 0.11 0.25 0.08 
Hg 2.223 12.43 + 0.88 0.99 + 0.37 0.61 0.21 
Hg 2.247 2.55 + 0.45 0.20 + 0.08 0.53 0.08 
Hə 2.406 11.90 + 0.96 1.43 + 0.61 NA NA 
Hə 2.413 20.60 2 13.70 2.47 + 1.94 NA NA 
Ho 2.423 10.76 + 1.10 1.29 + 0.56 NA NA 
Ho 2.437 7.12 + 0.90 0.85 + 0.37 NA NA 
Hə 2.454 5.86 + 0.84 0.70 a: 0.31 NA NA 
Hə 2.475 9.73 + 0.93 1.17 2 0.50 NA NA 
H2 2.499 4.61 + 0.79 0.55 + 0.25 NA NA 
Hə 2.550 4.54 + 0.84 0.54 + 0.25 NA NA 
H» 2.558 4.15 Æ 1.19 0.50 + 0.25 NA NA 
Hə 2.559 2.82 3E: 1:14 0.34 + 0.20 NA NA 


H, emission line fluxes and ratios with respect to H, 2.122 um (with flux (8.34+3.50)x10° Jy). Predicted ratios from the shock and fluorescence models are also given. 
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